Exact quantum spin liquids with Fermi surfaces in spin-| models 
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An emergent Fermi surface in a Mott insulator, an exotic quantum spin liquid state, was suggested 
by Anderson in 1987. After a quick support for its existence in spin-i Heisenberg model in a square 
lattice in a RVB mean field theory, pseudo Fermi surface was found only recently in an exactly 
solvable spin-| model by Yao, Zhang and Kivelson. We show that a minimal spin-| Kitaev model 
on a decorated square lattice exhibits a Fermi surface. Volume and shape of the Fermi surface 
change with exchange couplings or on addition of a 3 spin interaction terms. 

PACS numbers: 



Fermi surface emerges as a remarkable consequence of 
Pauli exclusion principle. It is a non-trivial organiza- 
tion of itinerant fermions in their ground states. So it 
has a natural place in describing a metal, where there 
are itinerant electrons. However, it was conjectured by 
Anderson[l|, that even in the absence of itinerant elec- 
trons a pseudo Fermi surface could arise from itinerant 
spins in a Mott insulator. Such a possibility was soon 
confirmed in a resonating valence bond (RVB) mean field 
theory for a spin-i Heisenberg antiferromagnct on a 
square lattice. In an exciting development, KitaevQ con- 
structed a spin-i model on a honeycomb lattice, where 
gapless noninteracting fermionic excitations emerge in a 
most natural fashion as an exact solution. There is no 
Fermi surface however, only a Fermi point. It is an im- 
portant question if one can change the model suitably 
and get a Fermi surface enclosing a finite volume of the 
Brillouin zone; and whether one can cause shape and 
volume change of the Fermi surface by tuning coupling 
parameters. 

On the experimental front, after years of efforts 
certain spin-i Mott insulating organic ET salts 0, and 
Na^^Os^, a 3 dimensional hyper kagome lattice anti- 
ferromagnet, exhibiting quantum spin liquid containing 
a pseudo Fermi surface have been synthesized. The spin 
Hamiltonian for the organic insulators, unlike the model 
studied in the present paper, has a SU(2) global spin ro- 
tational symmetry. The present model, which does not 
have a global SU(2) symmetry, has however features that 
should teach us about these real systems and systems 
describable by Kugel and Khomskii quantum compass 
mo del [3]. 

We show exactly that spin-^ Kitaev model defined on 
a suitably decorated square lattice (Fisher lattice) pos- 
sesses the desired psuedo Fermi surface in certain flux 
sector. We consider the translationally invariant model 
containing three types of couplings J x , 3 y and J z . We 
get a well defined fermi surface at the isotropic point: 3 X 
— 3 y = 3 Z = J. As we move away from this point the 
shape and volume of the Fermi surface change and even- 
tually disappear. A 3 spin interaction term also induces 
change in shape and volume of the Fermi surface, leading 
to an eventual disappearance of the Fermi surface. In a 
recent work Yao, Zhang and Kivelson (YZK) Q have sug- 
gested a spin- 1 model on a square lattice where they find 





FIG. 1: (a) decorated square lattice (Fisher lattice) with zero- 
flux in each plaquette. The gauge variable Uij (see text) is 
+ 1 if the arrow in the bond < ij > points from i to j. (b) 
dotted lines represent 3-spin interactions and the arrows on 
them represent the values of the variable Zij (see text). 



pseudo Fermi surface for a range of parameters. Nussi- 
nov and Ortiz |9j also discussed a class of models, whose 
spin representations include a spin-i system on a square 
lattice bilayer, with Fermi surface. These models may be 
viewed as non-trivial generalisation of the Kitaev model. 
The merit of our work is that even the simpler spin-i Ki- 
taev model, defined on a decorated square lattice(Fig. [T]) 
leads to a pseudo Fermi surface. The advantage of our 
construction is that it allows for possibility of Fermi sur- 
face even in certain three dimensional lattices [101 ] . Ki- 
taev model on a variety of interesting lattices was studied 
by [H. 

Ordered quantum antiferromagnetic states and their 
elementary excitations are well described by standard 
spin wave theories. However, these semiclassical theories 
are not capable of describing a quantum spin liquid state 
and their elementary excitations. This is because there is 
no manifest long range order describable by standard or- 
der parameters. It is here, RVB mean field theory broke 
new grounds in 1987 and introduced methods that were 
novel and unconventional. It is also interesting to note 
that RVB mean field approach, though yields approxi- 
mate results for isotropic Heisenberg antiferromagnet in 
2D, reproduces exact results for the Kitaev model 12j, an 
anisotropic spin model. 

In view of this we briefly review RVB mean field the- 
ory of pseudo Fermi surface of spin-^ Heisenberg anti- 
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ferromagnet on a square lattice. This theory differs from 
Anderson-Fazekas' [13| variational approch and was more 
ambitious: it developed a new formalism to study quan- 
tum spin liquid states and their excitations. It involved 
a key step of enlarging the Hilbert space by writing the 
spin operators in terms of electron operators constitut- 
ing the spin half magnetic moment and introducing bond 
singlet operators. In this representation the Heisenberg 
Hamiltonian takes the form, 

using the relation = J2 a b AoJap c iBi where r Q are 
the Pauli spin matrices and c's are the electron operta- 
tors that constitute the physical magnetic moment. Fur- 
ther b\- = ^5( c i| c J| — c Ii c ]t) ^ s t ^ ie hond singlet opera- 
tor. In the electron representation we have an enlarged 
Hilbert space of dimensions 4 . The local constraint, 
rij-f + riii = 1, reduces 4 N to 2^, the dimension of the 
physical Hilbert space. The Hamiltonian enjoys a local 
U(l) symmetry because the local electron number op- 
erator Hi = nj| + riii commutes with the Hamiltonian: 
[H,nj] = 0. Thus sectors containing same total number 
of particles N, but with N — 2M singly occupied sites, M 
doubly occupied and M unoccupied sites (M = 0, 1, 2, ... 
^r)are superselected. Unlike the Kitaev Model, enlarg- 
ing the Hilbert space does not allow for exact solvabil- 
ity of the isotropic Heisenberg antiferromagnet. However 
it suggests important approximation methods and solu- 
tions. 

A Bogoliubov-Hartree-Fock factorization leads to a 
mean field Hamiltonian: 

H mF = -JX;«^)6« + (hMj) ( 2 ) 

<«) 

and simple mean field solutions (6y) = Ao- After diago- 
nalising the mean field Hamiltonian we get the quasipar- 
ticle Hamiltonian: 

H mF ~ J I (cos k x + cos k y )\ a\ a a ka (3) 

The quasi particle energy vanishes on a square connect- 
ing points (±7r, 0), (0, ±7r) in the BZ. This defines the 
pseudo Fermi surface for spinon excitations. Unlike Ki- 
taev model, the RVB mean field analysis of the Heisen- 
berg Hamiltonian did not lead to 2 N gauge copies with 
identical energy spectrum. So one needs to Gutzwiller 
project the mean field ground state to the physical sub- 
space containing only single occupancies and then calcu- 
late the ground state and low energy properties. 

Now we study Kitaev model on the decorated square 
lattice (Fig. [I}. Kitaev model Hamiltonian is: 

H=-j x Y J ^ x j -j v Y J ^-JzY,^ (4) 

(*i>* {ij)v 



where (ij) x , (ij)y, (ij) z are the x, y and z-type nearest 
neighbor bonds shown in Fig. [T] There are two types of 
elementary plaquettes (square and octagon) in the deco- 
rated square lattice (Fig.[T]). There are two types of local 
conserved operators, S p — o\a 2 a^a\ for the squares and 
O p = o\a\G\a\a x a x a\a\ ii for the octagons. The pla- 
quette operators S p and O p commute among themselves 
and with Kitaev Hamiltonian. 

To solve for the spectrum of the Hamiltonian, we fol- 
low Kitaev and introduce 4 types of Majorana fermions 
on each sites, c a , a — 0,x,y,z, defined by the anti- 
commutation relations {c a ,c^} — 2S a g. Four Majorana 
(real) fermions make two complex fermions, making the 
Hilbert space 4 dimensional. This is similar to enlarge- 
ment of Hilbert space in RVB mean field theory. In fact, 
c a, s could be combined to form two complex fermions; 
c T = (c° + ic x )/2, c\ = (c° - ic x )/2, c L = (c* + ic z )/2, 

c\ = (c y -i(?)/2. 

Total number of spins in our system is 4N, where N 
is the number of unit cells. The dimension of Hilbert 
space of 4N spins is 2 4N . The enlarged Hilbert space has 
a dimension 4 4W = (\/2 x s/2 x V2 x V2) 4N . Hereafter 
we denote c° by Cj. State vectors of the physical Hilbert 
space satisfy the condition, 

A|*)phys = |*)phys, A = Ci Cfcfcf (5) 

The spin operators can then be represented by, erf = icicf 
where a £ {x,y,z}. When projected into the physical 
Hilbert space, the operators defined above satisfy the al- 
gebra of spin 1/2 operators, = 2ie a t, c a c 5ij . The 
Hamiltonian written in terms of the Majorana fermions 
is, 

H = ~ Y J «E iC ^(v)a C j' ( 6 ) 

a=x,y,z (ij) a 

with U(ij) a = icfCj. Following Kitaev we find that 
[H, M(ij)J = and it(y) a become constants of motion 
with eigen- values = ±1. The variables are 

identified with static (Ising) Z 2 gauge fields on the bonds. 
As Uuj) a = —u/ji) a we follow a definite ordering of the 
indices ij in specifying the value of U{ij) a , as shown in 
Fig. [TJ Kitaev Hamiltonian (Eq. [6]) has a local Z 2 gauge 
invariance in the extended Hilbert space. For practical 
purposes, the local Z 2 gauge transformation amounts to 
u (ij)a — > T i u (ij)a T ji with Tj = ±1. Eq.[5]is the Gauss law 
and the physical subspace is the gauge invariant sector. 

Eq. O with conserved uuj) a is the Hamiltonian of free 
Majorana fermions in the background of frozen Z 2 vor- 
tices or 7r-fluxes. Since Z 2 gauge fields have no dynamics, 
all eigenstates can be written as products of a state in 
the 2a( 4Ar ) dimensional Fock space of the Ci Majorana 
fermions and the (2)a( 4Ar ) dimensional space of Z 2 link 
variables. We will refer to the former as matter sector and 
the latter as gauge field sector. Gauge copies (eigen-states 
with same energy eigen-values) spanning corresponding 
extended Hilbert space are obtained by local gauge trans- 
formation on uujy 
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In the gauge field sector we have gauge invariant Z<i 
vortex charges ±1 (0 and 7r-fluxes), defined as product 
of uuj\ a around each elementary square and octagonal 
plaquette. These Z 2 fluxes (0, it) correspond to eigen 
values (-1, +1) of the plaquette operators S p and O p . 

For a given gauge flux sector, the dimension of the 
Fock space of the matter sector is 2 2N rather than 2 4JV . 
This is because we have a hopping Hamiltonian of Ma- 
jorana fermions rather than complex fermions which has 
27V complex fermion oscillators with energy eigen values 
> within the square Brillouin zone. The size of the BZ 
is where a is the periodicity of the lattice. 

We have analysed the spectrum of this model in dif- 
ferent flux sectors looking for pseudo Fermi surface. We 
succeeded in finding a Fermi surface in the sector, where 
each S p and O p has eigen value —1. This corresponds to 
Z 2 charge of +1 (zero flux) in all elementary plaquettes 
in the gauge field sector. By use of Lieb's theorem we 
find that the zero flux sector (or equivalently all S p and 
Op = —1) is not the minimum energy sector. However, 
we can make it a minimum energy sector by an addi- 
tion of a chemical potential term (S p + O p ), with 
j » 1, which commutes with the Hamiltonian. 

The gauge is chosen such that u/ij) a = +1 if, in 
Fig. []Ja), the arrow in the bond (ij) a point from i to 
j. The resulting lattice is periodic with each diamond 
forming the unit cell which is repeated along primitive 
vectors ni = (1,0) and n 2 = (0, 1). Using a new nota- 
tion Ci\ for the Majorana fermions, where i correspond 
to the unit cell and A the position within it, the Fourier 
transformation a \ (q) = J2i e~ tcl ' ri Cix/V2N transforms 
the Hamiltonian to 



H = 2 a x(l) A xA<lW(l)> 
1 

with {a\(q), at(fc)} = S q: kS\,^, a\(q) = a\(-q) and 
iA(q) = 2 



i JyO 



j x( jy 

iJ z a' 



—iJyCT 1 



(7) 



(8) 



where a=diag(e 192 , — e~ tqi ) and qi = q • n,. The skew- 
Hermitian matrix A(q) of size 2n x 2n with n = 2 has 
the property A T (q) = —A(—q), If Vj(g) is an eigenvec- 
tor of iA(q) with positive eigenvalues €i(q), then v*(— q) 
is also an eigenvector of iA(q) with eigenvalue — q)- 
Therefore, we could seperate the spectrum of iA(q) into 
m(q) positive eigenvalues and 2n — m{q) negative eigen- 
values with m(— q) = 2n — m{q). Then, the transfor- 
mation, di(q) = v*Aq)a {q), d' k (q) = v kj (-q)a j (q), i € 
{1, ..,m(q)}, k e {1, .., 2n — m(q)} diagonalises H, and 
the new operators are related by d\(q) — d'^—q). The 
Hamiltonian becomes 



H = 



-, m (i) 

q i—1 
m(q) 

J2T, e ^ k(?)di(?)-l/2 



(q) [4 (?)*(«) - 4(-q)4(-Q)] (9) 



(10) 




FIG. 2: Spectrum of iA(q) forJ x = J y = J z = 1 along 
(91,92) = (x,—x) (left) and along (91,92) = (as, as) . 
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FIG. 3: dependence of the size and shape of the Fermi surface 
(in a suitable gauge) on the coupling strengths J a when all 
plaquettes have zero flux. We fix J x — J y = 1 and J z takes 
values 0.5 and 1.9 respectively for the left and right figures. 



with the new operators obeying {di(q), dj(q')} = 5ij5 qq >. 

The ground state energy is given by — i J2 q 53i=i e i(<l)- 
For gapless excitations, at least one member of the spec- 
trum should be zero, the condition for which is given 
as Det(iA) = 0. It is readily seen that this is possible 
if the condition (J x — J z ) 2 < J 2 < (J x + Jz) 2 , which 



becomes J x < J v 



J z 1 Jy — Jx 



J z , J z — Jx "t" Jy 



when all J a are positive, is satisfied. Fig. [2] shows the 
eigenvalue spectrum €i(q) of iA(q) in Eq. [8] along lines in 
the (q x ,q y ) — (as, —as) and (q x ,q y ) — (x,x). The points 
= belong to the Fermi surface. One could see that 
these points are not symmetric about (0, 0) in this gauge. 
To exhibit the variation of size and shape of the Fermi 
surface with the parameters J a , we use a gauge which 
has the symmetry of the original lattice and the result 
is reported in Fig. [3l The volume of the Fermi surface 
becomes zero as we move out of the triangle defining the 
above mentioned inequality. 

To study topological properties of the gapless phase, 
characterised by Chern numbers, we added 3-spin inter- 
action terms (which is equivalent to adding an external 
magnetic field). Unlike the Kitaev model on a honey- 
comb lattice, this does not open a gap. Any non-trivial 
topological property of the pseudo Fermi sea remains hid- 
den at the moment. The 3-spin interaction, for the spins 
in the central square plaquette given in Fig. [TJa) can be 
written as K (erf erf erf + er|er|erf + a\a\a\ + a^ala^)- 
In terms of the Majorana fermions, the additional term 
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in the Hamiltonian becomes 



ff(3) 



iK ^ 

jk 



j kGjCki 



(11) 



where the sum is over second-nearest-neighbors and 
Zjk = —Zkj depends on the chosen guage u.y. In fig.HJb) 



these interactions are shown as dashed lines and Z«. 
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on the link (ij) if the arrow points from i to j. For the 
second-nearest-neighbors in the square plaquette there 
are two different 3-spin interactions that involves the 
same end spins. It turns out that they cancel each other 
in the gauge we have chosen. For example, in Fig. Q^a), 
spins 1 and 3 interact through crfo-f of and o% > o%o\ which 
in the chosen gauge adds upto zero when represented in 
terms of Majorana fermions. The effect of the 3-spin 
interactions is to modify A(q) by adding to it the matrix 



A^\q) = 2K 



Ai(g) A 2 (g) 

-AJ(«) A 2 fa) 

-A|fa) AJfa) 

-A£fa) -Aifa) 



,(12) 



where Aifa) = e iqi + e iq2 and A 2 fa) = er iqi + e iq2 . Fermi 
surface volume decreases with increasing K and vanishes 
asymptotically. 

Thus we have a Fermi surface as an exact solution of a 
non-trivial spin-i model in 2 dimensions. An important 
issue is the meaning of the area enclosed by the Fermi 
surface. Is there a Luttinger theorem that controls the 
area enclosed by the Fermi surface. Particularly surpris- 
ing is the fact that the volume of the FS changes with 
change in ratio of the coupling constants (Fig. [3]). 

We rationaise our result by comparing the band struc- 
ture of our problem with the problem of a complex 
fcrmion hopping on the same lattice in the same flux 
sector, the hopping amplitude between nearest neigh- 
bors being J a , namely we consider the Hamiltonian 

H = — J2 a JaJ2(ij) C l c j + ^- C - 1 R t ms case we have 
four single particle bands, because there are four sites 
per unit cell. Dimension of the Fock space is 2 4A \ The 
ground state corresponds to filling the bands upto zero 



energy, because of particle-hole symmetry. We have two 
hole like and two electron like bands. Two of them 
cross zero energy giving us two Fermi surfaces, one hole 
like and other electron like. Since we have two elec- 
trons (or two holes) per site in the state, Luttinger the- 
orem in this case efffectively demands only equality of 
the fermi area of the hole and electron pockets. Fur- 
ther, for the characteristic polynomials Pcf{qi, <? 2 , A) of 
the complex fermion hopping (Hermitian matrix) prob- 
lem and Pmf(qi, <]2, A) of Majorana fermion hopping 
(skew-Hermitian matrix) problem, in the gauge we have 
chosen earlier, are connected by the following relation 
Pcf fai, 92, A) = P M f fai + f , q_i + ir, A). There is also 
a two to one correspondence between the positive en- 
ergy eigen values of particle and hole excitation branches 
of a regular fermi sea and the positive energy complex 
fermion excitation branch of the Majorana fermi sea. 
That is, in the complex fermion problem we have de- 
generate positive energy excitations which are hole and 
electron branches. In the case of Majorana fermion they 
get identified, meaning that a complex fermion excitation 
is its own antiparticle. 

The present problem can also be solved using Jordan 
Wigner transformation, without enlarging the Hilbert 
space [lij • In this case the Jordan Wigner fermions have 
p-wave nearest neighbor pairing term with certain fixed 
pairing amplitude. What is interesting is that this pair- 
ing does not manage to open even partial gap in the 
Fermi surface. This is reminescent of the RVB mean 
field solution reviewed earlier, where an extended-S pair- 
ing leads to a Fermi surface. Further, nesting instability 
is imminent, as we have two Fermi surfaces with identi- 
cal shapes. This is also very similar to RVB mean field 
result for square lattice Heisenberg antiferromagnet. Ad- 
ditional perturbation will lead to nesting instabilities. It 
will be interesting to study this along the lines of work 
by Mandal, Sengupta and Shankar [15|. Equal time spin 
spin conrrelation function in this model also continues to 
vanish beyond nearest neighbors, as expected in all Ki- 
taev models [l6j. So what we have is a non trivial spin 
liquid with ultra short range spin-spin correlations that 
supports gapless Fermi surface excitations. 
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